Quantum Separability of the vacuum for Scalar Fields with a Boundary 
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Using the Green's function approach we investigate separability of the vacuum state of a massless 
scalar field with a single Dirichlet boundary. Separability is demonstrated using the positive partial 
transpose criterion for effective two-mode Gaussian states of collective operators. In contrast to the 
vacuum energy, entanglement of the vacuum is not modified by the presence of the boundary. 
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Recently theories and experiments in quantum in- 
formation science 0, 0] have been progressed signifi- 
cantly. This has led to interest in studying entan- 
glement in many-particle systems such as Bose Ein- 
stein condensations Q, Fermion systems 0, 0, lattice 
systems|j|, thermal boson systems 0,13 and even in the 
vacuum^ 0, [n|. Entanglement is now treated as a 
physical quantity, as well as a resource for quantum in- 
formation processing. On the other hand, scalar parti- 
cles such as the Higgs particles in the standard model 
are essential ingredients of modern physics. Although it 
is well known that the vacuums in some quantum field 
models have quantum nonlocalitv|l2j|. the existence of 
entanglement of the general vacuum is still controver- 
sial. From the particle physics viewpoint, everything 
in the universe is made of quantum fields and mass- 
less non-interacting quantum scalar field is the simplest 
and, hence, the most basic quantum field. Furthermore, 
scalar fields can be treated as an infinite-mode exten- 
sion of harmonic oscillators|l3j. Therefore, studying en- 
tanglement of the vacuum state of scalar fields is very 
important. For decades the nature of the vacuum, con- 
cerning the Casimir energy ^IE]> in a bounded space 
has been extensively studied using the Green's function 
method[l^|. In this paper, we use this Green's function 
approach to investigate the quantum entanglement of the 
vacuum for non-interacting scalar fields with an infinite 
Dirichlet boundary (i.e., fields are constrained to vanish 
at the boundary). Since the scalar fields are continuous 
variables we need a separability criterion for continuous 
variables. Recently there has been a renewal of interest 
in entanglement of Gaussian states in the context of con- 
tinuous variable quantum information, and entanglement 
measures such as the puritv|l7| or the negativity has been 
suggested^ f° r two-mode Gaussian states. In this work 
we shall average the fields over two tiny boxes to make 
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the vacuum state an effective two-mode Gaussian state. 

We start by introducing the massless (real) Klein- 
Gordon fields, which are described by a Lagrangian den- 
sity, 



V, (a = t,x,y,z), 



(1) 



which leads to the equation of motion d a d a <fi = in 
unbounded four-dimensional Minkowski spacetime. In 
this free-space (without a boundary) the quantum field 
operator of the scalar field can be expanded as 



(j>(x, t) = 




a(k)e 



ik-x—itJkt 



+ H.C. 



(2) 



where Uk = \k\. The equal time commutation relations 
of the scalar field are 



(x,t),(j>(x,t)] = 0,[Tr(x,t),ir(x',t)}=0, 
(x, t), tt(x' , t)] — iS(x — x'), 



(3) 



where n(x,t) = dt<p(x,t) is the momentum operator for 
4>{x,t). Consider a vector of the field and its momen- 
tum operator at two points x and x' at a time t, i. e., 
£ = {4>(x, t), 7r(x, t), 4>(x' , t), 7r(x', t)). The vacuum for the 
scalar field (|0)) has a variance matrix having the follow- 
ing form: 

V af3 = ^(Q\{A£ a ,&&}\0), (4) 

where {A, B} = AB + BA and A£ Q = £ Q - (0|£ Q |0) with 
(0|£ Q |0) = in this paper. This matrix represents zero- 
temperature quantum fluctuation of the vacuum. Us- 
ing the equal time commutation relations one can easily 
note that (0|&6|0) = (0\<f>(x, i)0(f, *)|0) = (016610), 
hence V13 = V31. For systems having time-translational 
symmetry, like ours, dt(0\4>(x,t)(f>(x,t)\0) = 0. There- 
fore, (0|£i&|0> = (O|0(i?,t)7r(f,t)|O> = -(0|66|0), and 
V12 = V21 = 0. Hence the variance matrix has a following 
form 



V 



a c 0- 

b d 

c a' 

d b' 



(5) 



2 



where a = (O|{0(z, t), (p(S, t)}\0)/2, b = 
(0|{7r(f,t),7r(j?,t)}|0)/2, a' = (0\{cf>(x>,t),<t>(x>,t)}\0)/2, 
and so on. Note that for the free-space field operator 
(Eq. J5J) a and a' term diverge. The divergence, 
however, can be made disappear if we consider a scalar 
field with a boundary. The effect of the boundary can be 
calculated by subtracting a free-space Green's function 
from the Green's function with a boundary (See Eq.JSJ 
and below). Divergent terms which are in the both 
Green's functions cancel each other. The components 
of V can be calculated from the Green's function called 
the Hadamard's elementary function 0] 



G(x,t;x' ,t') = (O|{0(:r,t),(Mx',i')}|O}. 



(6) 



Once we know the Green's function we can calculate the 
components of the variance matrix V from it, i.e., 



c(x, x') 



d(x, x") 



1 



(0|{<Mx,t), </>(£', i)} 1 0} 



(7) 



lim ~2&{x, t; x', i'), 

±(o\{M(x,t),d t 0(x',t)}\o) 

lim d t d t ' -jG(x, t; x' , t'). 



Then, a = lirng/^ c(x 7 x'), b = lim^/^j d(x,x') and so 
on. 

With the boundary the mode expansion in Eq. © is 
no longer valid. Instead one can directly work with a 
modified Green's function to obtain V as follows. First 
let us choose two points x = (x, y, z) and x' = (x' ', y' , z'). 
Using the method of images [2(| one can obtain a scalar 
field Green's function Gb with a Dirichlet boundary po- 
sitioned at z = 



Gb — Gn — 



1 



!vr 2 (-(t 



t'f +r 2 + {z + z'f 



(8) 



where the free-space Green's function is 
G = [2ir 2 (-{t-t'f + r 2 + {z-z'f^j]- 1 and 

r 2 = (x — x 1 ) 2 + (y — y') 2 . By subtracting the free-space 
Green's function from Gb we obtain a regularized 
Green's function; G = G b — Gq which is not divergent. 
From G one can obtain the components of V using Eq. 
0; 



{a,b,a',b',c,d) = r^o(r-4, 



-1 



1 



2tt 2 '8 Z 2 ' 16 z 4 ' 8z> 2 ' 16 z' 4 ' 
=1 1 ) 

2 >■ 



(9) 



r 2 + (z + z'f) (r 2 + {z + z') 2 )' 



Now we discuss separability of the vacuum. Since 'if and 
only if ' separability test for infinite-mode states are un- 
known, we need to reduce the infinite-mode states to 
effective two-mode Gaussian states. We follow the ap- 
proach in ref. 13], i.e., we spatially average the field 



B' 



B 



FIG. 1: We average fields over two tiny boxes B and B' 
centered at x and x" , respectively. The Dirichlet boundary is 
at z = 0. 



operator over two tiny boxes centered at x and x' , respec- 
tively (See Fig. 1). Defining these collective operators is 
reasonable, because in a physical situation probes always 
have finite spatial resolution. We, however, do not av- 
erage but integrate the momentum operators within the 
box, since momentum is additive. Thus the collective 
operators are 



I = ($(x,t),n(f,<),$(f,<),n(f ,*)) 

= (^3 / d 3 y(f>{x + y,t), f d 3 yjT(x + y,t), 



(10) 



£3 / d^y'Mx' 



y',t), I d 3 y'Tr(x" + y>,tj). 

' B' 



Here J B d 3 yf(x + y) denotes an integration of f(x) 
over a box B centered at x with volume L 3 . Then, the 
commutation relations in Eq. reduce to the canonical 
commutation relations for the collective operators 0] 



a/3, 



where 



",/ 


0" 


,J = 


' l" 





J 


-1 



(11) 



(12) 



Similarly we define an variance matrix for the collective 
operators V a p = ^(0\{^ a , 1/3 }|0). 

Now we show that, using the fundamental theorem of 
calculus [2l|. one can calculate the variance matrix V for 
the two-mode Gaussian states with those of V in the limit 
L — ► (but still non-zero). For example, 

^13 = lim i(0|{$(f,t),$(^',i)}|0) 

= lim ^ J d 3 y j dV(O|{0(f + y, t), cj>{x' + ff t)}\0) 

= ^(0\{4>(x,t),<t>(x>,t)}\0) = V 13 = c. 

This is possible since the integrand of the second line (i. 
e., c term in Eq. ©) is continuous on all the region 



(13) 



3 



z > 0. ( Due to the boundary, it is enough to consider 
only the half-space z > 0.) Similarly 



% A = limi 6 - J d 3 yj^d 3 y'(0\Mx + y ) t)(U) 



tt( z' + ^, t)}|0) = L 6 V 2A = L 6 d. 



Hence, 



V 



a c • 

L 6 b L 6 d 

c a' 

L 6 rf L 6 b' 



A G 
G T B 



(15) 



Obviously, our formalism does not work for L = 0, but for 
L infinitesimally small the above equation gives an exact 
result. The separability criterion we use in this paper 
is positive partial transpose (PPT) criterion]^, H3| for 
two- mode Gaussian states |24j which is equivalent to 



F = E - (- + AdetV) < 0, 



(16) 



where E = detA + detB — 2detG. By inserting the com- 
ponents into Eq. (|16() we obtain 



1 

L 6 

2 2 7T 4 

L 12 



(17) 



1 



1 



1 



2 7 z 6 
1 



2 7 z /6 [ r 2 + (z + z') 2 ] 3 
1 



2 6 tt 8 L 2 12 z 6 z' 6 2 4 zV 2 (r 2 + (z + z') 2 ) 4 
1 1 



< 



2 8 zV 4 (r 2 + (z + z') 2 ) 2 (r 2 + (z + z') 2 ) 6J 
1 

~~4' 



One can obtain this inequality using two absolute in- 
equalities, 1/(X + Y) n < (1/X n + l/Y n )/2 n+1 and 
X 3 + Y 3 > XY 2 + YX 2 (X, Y > 0) on the two square 



brackets, respectively. The maximum value is achieved 
when x = x' (r = 0). Therefore, the effective two- 
mode scalar field vacuum with a boundary is PPT and, 
hence, separable, when it is described by the variance 
matrix V and the regularized Green's function G. Our 
results, however, do not rule out the possibility of entan- 
glement of the vacuum for free-space scalar fields, since 
we have subtracted the free-space Green's function which 
might give entanglement. What our results really im- 
ply is that the presence of a single boundary does not 
change the separability or the entanglement which the 
vacuum for the free-space scalar fields may have, while 
the presence of the boundary modifies the vacuum energy 
to -l/167T 2 z 4 [I3. 

Entanglement of the scalar field could be experimen- 
tally tested by the scheme with trapped ionsQ or Bose- 
Einstein condensates |25j| . Our approach provides a new 
method using the Hadamard's elementary function to 
investigate entanglement of the quantum fields vacuum 
within a bounded space. It will be interesting to inves- 
tigate how the number of boundaries and properties of 
the fields such as interactions, masses, charges, and spins 
change the results. 



After completion of our work we found that there ap- 
pears a paper about spatial entanglement of free thermal 
bosonic fields (quant-ph/0607069). 
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